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The electron-transport properties of adatom-graphene system are investigated for different spatial 
configurations of adsorbed atoms: when they are randomly-, correlatively-, or orderly-distributed 
over different types of high symmetry sites with various adsorption heights. Potassium adatoms 
in monolayer graphene are modeled by the scattering potential adapted from the independent self- 
consistent ab initio calculations. The results are obtained numerically using the quantum-mechanical 
Kubo-Greenwood formalism. A band gap may be opened only if ordered adatoms act as substi¬ 
tutional atoms, while there is no band gap opening for adatoms acting as interstitial atoms. The 
type of adsorption sites strongly affect the conductivity for random and correlated adatoms, but 
practically does not change the conductivity when they form ordered superstructures with equal 
periods. Depending on electron density and type of adsorption sites, the conductivity for corre¬ 
lated and ordered adatoms is found to be enhanced in dozens of times as compared to the cases 
of their random positions. These the correlation and ordering effects manifest weaker or stronger 
depending on whether adatoms act as substitutional or interstitial atoms. The conductivity ap¬ 
proximately linearly scales with adsorption height of random or correlated adatoms, but remains 
practically unchanged with adequate varying of elevation of ordered adatoms. Correlations between 
electron transport properties and heterogeneous electron transfer kinetics through potassium-doped 
graphene and electrolyte interface are investigated as well. The ferri-/ferrocyanide redox couple is 
used as an electrochemical benchmark system. Potassium adsorption of graphene electrode results 
to only slight suppress of the heterogeneous standard rate constant. Band gap, opening for ordered 
and strongly short-range scatterers, has a strong impact on the dependence of the electrode reaction 
rate as a function of electrode potential. 

PACS numbers: 72.80.Vp, 81.05.ue, 82.20.Pm 


I. INTRODUCTION electric-dipole moment, magnetic moments of an isolated 

atom and total graphene-adatom system, etc > — 


Adsorbed atoms and molecules are probably the most 
important examples of point defects in the physics of 
graphene A In addition to remarkable intrinsic electronic 
and mechanical properties of pure graphene, its struc¬ 
ture and properties can also be modified and controlled 
by adsorption and doping of atoms and molecules. That 
is why last few years studies of atom adsorption of 
both metallic*^— and nonmetallic42r— adsorbates on 
graphene attract a considerable attention. Overwhelm¬ 
ing majority of theoretical and computational studies 
of adatom-graphene systems deal with first-principles 
density-functional calculations, which require high com¬ 
putational capabilities, therefore the size of graphene 
computational domains in these calculations are mostly 
limited to periodic supercells and localized fragments 
containing a relatively small number of atoms (sites). 
Nevertheless, the first-principle study is suitable and 
fruitful, and therefore prevalent now, for calculation 
of energetic, structural, and magnetic parameters: ad¬ 
sorption (binding) energy and height of adatoms, dif¬ 
fusion (migration) barrier energy, in-plane and vertical 
graphene-lattice distortion amplitude, charge transfer, 


Because of the hexagonal symmetry of the graphene 
lattice, possible adsorption sites for a single atom can 
be reduced into three types with high-symmetry posi¬ 
tions: so-called hollow center (H- type), bridge center ( B- 
type), and (a)top (T-type) adsorption sites as illustrated 
in Fig. [T] The most favorable (stable) adsorption site is 
determined by placing the adatom onto these three ad¬ 
sorption sites, and each time by optimizing structures to 
obtain minimum energy and atomic forces; as a result, 
the highest binding (adsorption) energy of adatom cor¬ 
responds to its the most favorable site. Analysis of the 
density-functional-theory-based studies^ - — covering al¬ 
most all the periodic table, yields: (i) for metals, the 
most stable adsorption sites are the H- sites, followed by 
the .B-sites, and then the T-sites, although the energy 
differences between the H and B or T sites are very 
small for the alkali and group-III metals, particularly for 
potassium (see Table (H), which we regard as an example 
of adsorbate in the present study; (ii) for both metals and 
nonmetals, adsorption heights for more favorable sites are 
lower as compared with heights of the lesser favorable ad¬ 
sorption sites. 
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Figure 1: (Color online) Typical configurations of adatom-graphene system: top (left) and perspective (right) views of graphene 
lattice with hollow center ( H ), bridge center ( B ), and (a)top (T) adsorption sites. 


Table I: Literature data on calculated adsorption energies and 
heights for K adatoms occupying hollow ( H ), bridge ( B ), and 
top ( T ) adsorption sites in graphene. 


Calculated parameter 

Adsorption 

site 

H- type 

B-type 

T-type 


0.785 a 

0.726 a 

0.720 a 

Adsorption energy [eV] 

0.802 b 

1.461 c 

0.739 b 

1.403 c 

0.733 b 

1.405 c 


0.810 d 




2.62 a 



Adsorption height [A] 

2.60 b 

2.52 c 

2.67 b 

2.59 c 

2.67 b 

2.55 c 


2.58 d 


a d References MM (respectively). 


Data of Table Q] for K adsorption on graphene read 
that values of adsorption energy reported in the litera¬ 
ture disagree by as much as almost 100%, while adsorp¬ 
tion heights differ by up to 5%— Similar inconsistencies 
of the literature data occur also for other periodic-table 
elements. For example, Cu and Sn prefer T-site bonding 
(at the heights of 2.12 A and 2.82 A, respectively) ac¬ 
cording to Refs. I 22 I and fd8l. while R-site bonding (2.03 A 
and 2.42 A) in accordance with Refs. [33j and [1(1 O 11 
the one hand, such discrepancies in determination of the 
energy stability of adsorption sites has resulted in a con¬ 
troversy and questions concerning the accuracy of theo¬ 
retical models (calculations) used in those studies. On 
the other hand, this motivates us to study how the po¬ 
sitioning of adatoms on each of H , R, and T site types 
affects the transport properties of graphene in compari¬ 
son with the cases of their location on two other types of 
the sites. 


Distributions of adatoms over the H, B , or T graphene- 
lattice adsorption sites are not always random, as it is 
usually in three-dimensional metals and alloys, where 
adatoms are introduced by alloying, which is generically 
a random process— Diluted atoms may have a tendency 
towards the spatial correlation^ or even ordering! 50 ! 54- — 
Moreover, since graphene is an open surface, (ad)atoms 
can be positioned onto it with the use of scanning 
tunneling^ or transmission electron^ microscopes al¬ 
lowing to design (ad) atomic configurations as well as or¬ 
dered (super)structures with atomic precision. Recently, 
several ordered configurations of hydrogen adatoms on 
graphene have been already directly observed by scan¬ 
ning tunneling microscopy in Ref. 15ft 

Though many properties of atom adsorption onto 
graphene have been extensively studied in many works, 
there is still no one paper on how such a variety of the 
spatial arrangements of adatoms (viz., their random, cor¬ 
related, and ordered distributions in the H, B, and T 
types of bonding with varying adsorption heights) influ¬ 
ences (if any) on electron transport in graphene. Such a 
problem formulation arises in context of the possibility 
to consider (ad) atomic spatial configurations as an ad¬ 
ditional tool for modification and controlling graphene’s 
transport properties. 

Another part of our paper deals with attempt to detect 
adatom-mediated correlation between electron transport 
and electrochemical properties of graphene. Under¬ 
standing of its electrochemical properties, especially the 
electron transfer kinetics of a redox reaction between 
graphene surface (electrode) and redox couple in elec¬ 
trolyte, is essential^ for its potential in energy conversion 
and storage to be realized , 62 ! 63 as well as opens up inter¬ 
esting opportunities for using graphene as an electrode 
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material for field effect transistor s 64 : 65 and electrochem¬ 
ical senorsj£& 67 To examine the heterogeneous electron 
transfer kinetics at highly oriented pyrolytic graphite 
(HOPG) and glassy carbon (GC) electrode, several elec- 
troactive species were usedi^ Results show that electron 
transfer is slower at the basal plane of HOPG than at 
the edge plane. The kinetics of the electron transfer is 
enhanced after electrode pretreatment. However, in epi¬ 
taxial graphene, only a part of the surface is electroactive, 
even after electrochemical pretreatmentExperimental 
results confirmed the belief that point and edge defects 
as well as oxygenated functional groups can mediate elec¬ 
tron transfer^ Contrary to the traditional view, high- 
resolution electrochemical imaging experiments have re¬ 
vealed that electron transfer occurs at both the basal 
planes of graphite as at the edge sites^ To examine these 
discrepancies, we calculate the electron transfer kinetics 
at graphene with randomly-, correlatively-, and orderly- 
adsorbed atoms described by scattering potential mani¬ 
festing both short- and long-range features, and also use 
strongly short-range scattering potential. Results show 
that electron transfer still occurs for adsorbed graphene. 

The rest of the paper is organized as follows. Section II 
consists of two subsections containing models for electron 
transport and transfer. In the first subsection, we for¬ 
mulate the Kubo-Greenwood-formalism-based numerical 
model for electron transport in graphene, which is appro¬ 
priate for realistic graphene sheets with millions of atoms. 
The size of our computational domain is up to 10 mil¬ 
lions of atoms that corresponds to ~ 500 x 500 nm 2 . The 
second subsection encloses the basic model we use to cal¬ 
culate the rate constant of electron transfer between solid 
(graphene) electrode and redox couple in electrolyte us¬ 
ing the Gerischer-Marcus approach. Section III presents 
and discusses the obtained results. Finally, the conclu¬ 
sions of our work are given in Sec. IV. 


II. MODELS 
A. Electron transport 

To investigate the charge transport in adatom- 
graphene system, an exact numerical technique within 
the Kubo-Greenwood formalism,' 2 ~— which captures all 
(ballistic, diffusive, and localization) transport regimes, 
is employed. Within the framework of this approach, the 
energy (E) and time (t) dependent diffusivity D(E , t)2£ is 
governed by the wave-packet propagation: 7 VSi D(E, t) = 
(A X 2 (E, t))/tr& where the mean quadratic spreading of 
the wave packet along the direction x reads asZ^r— 


(A X 2 (E,t)) 


Tr[(Af(t) — X(0)) 2 S(E — H)\ 
Ti[8(E - H)} 


( 1 ) 


with X(t) = U'(t)XU(t )—the position operator in the 
Heisenberg representation, U(t) = e ~ lHt ^ h —the time- 


evolution operator, and a standard p-orbital nearest- 
neighbor tight-binding Hamiltonian H i o 94 i 95 


H = + EiVicJa, ( 2 ) 

where cj (cj) is a standard creation (annihilation) opera¬ 
tor acting on a quasiparticle at the site i. The summation 
over i runs the entire honeycomb lattice, while i' is re¬ 
stricted to the sites next to i; u = 2.7 eV is the hopping 
integral for the neighboring C atoms occupying i and i' 
sites at a distance a = 0.142 nm between them; and V) 
is the on-site potential defining scattering strength on a 
given graphene-lattice site i due to the presence of impu¬ 
rity adatoms. The impurity scattering potential plays a 
crucial role in the transport model we use at hand. 

For adatoms located on H-type sites (see Fig. [TJ, the 
impurity scattering potential in the Hamiltonian matrix 
is introduced as on-site energies V t varying with dis¬ 
tance r to the center of a hexagon on which the adatom 
projects according to the potential profile V = V(r) < 0 
in Fig. [2ja) adapted from the self-consistent ab initio 
calculations^ for K adatoms on the height h = 2.4 
A over the graphene surface. As the fitting22 shows, 
this potential is far from the Coulomb- or Gaussian-like 
shapes commonly used in the literature for charged impu¬ 
rities in graphene, while two-exponential fitting exactly 
reproduces the potential. Such a scattering potential 
presents both short- and long-range features^ 9 although 
its short-range characteristics become rather stronger for 
adatoms that are nonrandom (correlated and ordered) in 
their spatial positions. Transforming scattering potential 
V = V(r) into its dependence on distance from the lattice 
site directly to adatom, V = V(l), where l = \/r 2 + h 2 
as demonstrably from Fig. [T] one can obtain its depen¬ 
dence on both r and h,V=V (r, h), which is presented in 
Fig. (21(b). As follows from Fig. (21(b), if r = a and h = 2.4 
A = 1.69a, V = —0.37 u, which agrees with Fig. (21(a). 

For adatoms positioned on B- and T-type sites 
(Fig. CD, we use the same scattering potential V = V(r) 
as in Fig. (21(a) with difference that r denotes distance 
from the lattice site to the middle of a C-C bond and to 
a C atom, respectively. Strictly speaking, V = V(r) for 
adatoms on H-, R-, and T-type sites (Fig. [[]) should be 
different, however just approach of the same scattering 
potential for these three types of adatom locations allows 
us to reveal manifestation of configurational effects in the 
transport properties of graphene we are interested in the 
present study. 

In case of correlation, adatoms are no longer consid¬ 
ered to be randomly located. To describe their spatial 
correlation, we adopt a model^^ using the pair distri¬ 
bution function p(R^ — Rj) = p(r ): 


p(r) 


0, r <r 0 
1, r > r 0 


(3) 


where r = |Ri — Rj| is a distance between the two 
adatoms, and a correlation length tq defines minimal dis- 
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Figure 2: (Color online) Scattering potential for potassium adatoms with (a) fixed adsorption height h = 2.4 A and (b) varying 
h. Here, ab initio calculations (•)££ are fitted2£ by different functions, viz. Gaussian (V = Ue~ r with fitting parameters 
U = —0.37u and £ = 2.21a defining a potential height and an effective potential radius, respectively), Coulomb (V = Q/r with 
Q = —0.36 ua), and two-exponential (V = Uie~ r ^ 1 + L^e -7 ^ 2 with Ui = —0.45u, £i = 1.47a, U 2 = —0.20u, £2 = 2.73a); r is 
a distance from the projection of adatom to the lattice site as shown in Fig. [T] 


tance that can separate any two of them. Note that for 
the randomly distributed adatoms, ro = 0. Although 
the correlation length ro is found to be insensitive to 
impurity (potassium) density*** the maximal correlation 
length ro max depends on both relative adatom concentra¬ 
tion and positions of adsorption sites as given in Table HT1 
In our calculations for uk = 3.125% of correlated (potas¬ 
sium) adatoms, we chose rn = r^' B =7 a for hollow and 
bridge sites, and = 5a for top sites. 

In case of adatom ordering, we consider superlattice 
structures in Fig. [3) where the relative content of ordered 
(potassium) adatoms is the same as for random and cor¬ 
relation cases, ni< = 1/32 = 3.125%. This structures 
form interstitial [Fig. (3|left) and Fig. [S]) center)] or sub¬ 
stitutional [Fig. Eli'ight)] superstructures, where distri¬ 
bution of adatoms over the honeycomb-lattice interstices 
or sites, respectively, can be described by the single-site 
occupation-probability functions derived via the static 
concentration wave methodA— ~ 104 In the computer im¬ 
plementation, riK = 3.125% of potassium adatoms oc¬ 
cupy sites within the same sublattice and can be de¬ 
scribed via a single-site function: 


P(R) = 


1, Til + 712 = 4Z 
0, otherwise 


(4) 


where 712 , and Z belong to the set of integers, n\ and 
n .2 denote coordinates of sites in an oblique coordinate 
system formed by the basis translation vectors ai and 
a 2 shown in Fig. [3] and R denotes origin position of the 
unit cell where the considered interstice [Fig. [deleft) and 
Fig. [3][center)] or site [Fig. fright)] resides. 

The dc conductivity a can be extracted from the dif- 
fusivity D(E,t), when it saturates reaching the maxi¬ 
mum value, linq^oo D(E, t) = D max (E), and the dif¬ 


Table II: Relation between the relative concentration of im¬ 
purity adatoms (rii) occupying H-, B-, or T-type sites (see 
Fig.m and the largest correlation distance (ro max ) expressed 
in units of the lattice parameter a = 0.142 nm. 


Site 

Ui 

0 . 5 % 

1% 2% 3% 4% 5% 

H, B 

Jrl r I 

r o^. x M 

18 

13 

9 7 6 5 

T 

r o ma .x M 

13 

9 

6 5 4 3 


fusive transport regime occurs. Then the semiclassical 
conductivity at a zero temperature is defined a s 81 ’ 82 


a = e 2 p(E)D max (E), (5) 

where —e < 0 denotes the electron charge and p(E) = 
p/Q = Tr[6(P — H)\/Sl is the density of sates (DOS) per 
unit area O (and per spin). The DOS is also used to 
calculate the electron density as n e (E) = p(E)dE — 
fhorMo where rii ons = 3.9 • 10 15 cm -2 is the density of 
the positive ions in the graphene lattice compensating 
the negative charge of the p-electrons (at the neutrality 
(Dirac) point of pristine graphene, n e {E) = 0). Com¬ 
bining the calculated n e (E) with <r{E), we compute the 
density dependence of the conductivity a = <j(n e ). 

Note that we do not go into details of numerical calcu¬ 
lations of DOS, D(E,t), and a since details of the com¬ 
putational methods, we utilize here (Chebyshcv method 
for solution of the time-dependent Schrodinger equation, 
calculation of the first diagonal element of the Green’s 
function using continued fraction technique and tridiag- 
onalization procedure of the Hamiltonian matrix, aver¬ 
aging over the realizations of impurity adatoms, sizes of 
initial wave packet and computational domain, boundary 
conditions, etc.) are given by Radchenko et alJ£ 
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Figure 3: (Color online) Top view of graphene lattice with ordered adatoms resided on hollow (left), bridge (centre), and top 
(right) sites. 


B. Electron transfer 


To calculate the heterogeneous rate constant of elec¬ 
tron transfer from the reduced form of redox cou¬ 
ple to graphene electrode, we used Gerischer-Marcus 
modelii££c^i In this model it is assumed that electron 
transfer between solid electrode and redox couple in elec¬ 
trolyte is much faster than reorientation of the solvent 
molecules (diabatic representation). As a result, the rate 
constant of the electrode reaction depends only on the 
electron DOS in the solid and the distribution energy 
levels of the reduced (oxidized) form, Wf. ed ( q x ) in the 
solution. If the vacuum energy as a reference energy 
level is chosen, the electrochemical potential of electrons 
occupying energy levels of ions, fL e ,ie dox, is equivalent 
to the Fermi level of the redox couple in the solution, 
Ep, redox' 112 As oxidized and reduced form interact with 
surrounding polar solvent in a different way, energy levels 
of oxidized and reduced form are shifted each other by 
2A, where A is the reorganization energy. 

In our calculations, we used the Gaussian distribution 
of the electronic states of the reduced form given byiSS 


P(E) 


1 

V4fc B TA 


exp 


(E Ep,red ox 

4k B TX 



( 6 ) 


where k B is the Boltzmann constant, T is the absolute 
temperature. The Fe(CN) 3 _ ^ 4 ~ redox couple has been 
chosen as a benchmark system. For ferri-/ferrocyanide 
redox couple, the A value ranges between 0.5 eV and 
1.0 eVpii^ In our calculations we have used intermediate 
value of 0.71eViiili 

Dependence of the cathodic reaction rate on the elec¬ 
trode potential k c (V ) is given by the integra l 109 ! 110 

k c oc J [1 - f{E , V)\ DOS (E, V) P{E) dE. (7) 

Here, DOS(i?, V) = BOS(E - eV), and f{E,V) = 
f(E — eV), where f(E) is the Fermi Dirac distribution. 
To determine the position of electron bands of graphene 
electrode in relation to the Gaussian distribution of en¬ 
ergy levels of the reduced form, vacuum energy has been 
chosen as a reference. The value of Ep in relation to the 


vacuum energy is equal to the work function, which has 
been determined experimentally for mono- and bilayer 
graphene using Kelvin probe force microscopy giving the 
Ep[vs. vacuum] = —4.6 eVJ^ 

For the Fe(CN) 3 _ ^ 4_ redox couple, we used the 
Ef, redox [ws. vacuum] value determined from the half wave 
potential obtained by cyclic voltammograms. Measure¬ 
ments carried out on epitaxial graphene and HOPG give 
value ranging from V) , 2 = —0.025 V vs. Ag/AgCl to 
V 1/2 = 0.268 V vs. Ag/AgCl^ According to Ref. Ill5l 
zero potential of the Ag/AgCl reference electrode is 
shifted in relation to the vacuum potential by —5.04 V. 
Assuming the half wave potential value to be 0.20 V, 
we found the value of -Ef, redoxes. vacuum] = —4.84 eV. 
Thus, we have assumed in our model that the Fermi level 
of the Fe(CN) 3 _ ^ 4_ redox couple is shifted in relation to 
the Fermi level of the graphene by —1.27 eV. 


III. RESULTS AND DISCUSSION 

As it was mentioned in Secs. I and II, we consider 
potassium as an example of adsorbate in the present 
study. The most energy favorable adsorption sites for K 
dopants in graphene are H-type sites as listed in Table |T] 
Therefore results obtained at potential in Fig. [U]a) and 
17-type sites are appropriate for K adsorbed graphene 
first of all. Results obtained for B- and T-type sites can 
be associated with K adsorbate in a model assumption 
for revealing manifestation of configurational effects in 
electron transport. 

Figure @Ja) shows the DOS and the electron density 
n e = n e (E ) for graphene with tik = 3.125% of random, 
correlated, and ordered potassium adatoms, which are 
described by the scattering potential in Fig. [5]]a) and 
are distributed over the H-type adsorption sites. DOS- 
curves for B- and T-type are similar to those shown 
in Fig. [4ja) with difference that Fermi level in case of 
T sites is shifted more far with respect to E = 0 to 
the (left) side of negative E < 0—energies of holes in 
our denotations. The Dirac (neutrality) point shifts to¬ 
wards negative energies (gate voltage) due to electron 
(n-type) doping dictated by the asymmetry (negativity) 
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Figure 4: (Color online) (a), (b) Density of states (DOS) for (a) potential in Fig. UK a) and for (b) short-range Gaussian 
potential V = Ue~ r ' 4 (with potential height U = —2 u and effective potential radius £ = 0.5a), where (a) 3.125% of random, 
correlated, or ordered adatoms occupy hollow sites, while (b) 3.125% (stoichiometry 1/32), 6.25% (1/16), and 12.5% (1/8) of 
ordered adatoms reside on top sites, (c), (d) Rate constant (fc c ) of the cathodic reaction of oxidation of Fe(CN)g~ at mono-layer 
graphene electrode for respective (a), (b) DOS. Insets in (a) and (b) show the same as in the main panels, but with another 
scales. As a reference, DOS (a) and rate constant (c), (d) for pure graphene electrode are shown. 


of the scattering potential. The calculated DOS-curves in 
Fig. 0](a) for random, correlated, and ordered K adatoms 
are similar with two differences take place in the case 
of ordering: (i) peaks (fluctuations), which appear close 
to E/u « —3 at correlation, begin to manifest them¬ 
selves in all energy interval (weakly away from the re¬ 
gions of the van Hove singularities and \E/u\ « 3, but 
stronger close to them); (ii) at a Fermi energy level, the 
DOS drops to zero (but even small band gap does not 
open). Appearing of the peaks (fluctuations) in DOS is 
due to the periodicity of the scattering-potential distri¬ 
bution describing ordered positions of adatoms on the 
sites of interstitial [Figs. 0](lcft) and [Sjcenter)] or sub¬ 
stitutional [Fig. Enright)] graphene-based superstructure. 
Additional calculations^ [see also inset in Fig. [41(b)] show 
that the peaks become stronger and even transform into 
discrete energy levels with broadening as impurity con¬ 
centration and/or periodic potential increase. 

Positioning of ordered adatoms on the T-type sites 
[Fig. Enright)] makes possible band gap opening, which is 
clearly seen in Fig. [4Kb) for a strongly short-range, e.g ., 
Gaussian potential with very small effective potential ra¬ 


dius £ < a. The band gap is induced by the periodic 
potential leading to the ordered distribution of adatoms 
directly above the C atoms belonging to the same sub¬ 
lattice, thus breaking of symmetry of two graphene sub¬ 
lattices. Note that adatoms on T sites act as substitu¬ 
tional point defects— -impurities or vacancies—which can 
also induce the band gap opening if they are distributed 
orderl y 54 : 89 ! 116- -^ or belong to the same sublattice even 
being randomly located i 119 : 120 However, we did not ob¬ 
serve the band gap appearing if ordered adatoms reside 
on H and B adsorption sites (thereby act as interstitial 
dopants) as it is reported by Cheianov et al. for adatoms 
occupying H— and B — sites. We attribute the lack of 
the band gap opening (when ordered adatoms occupy H 
and B sites) to the absence of the breaking of global lat¬ 
tice symmetry in these cases. 

Obtained densities of electronic states enter into 
Eq. 0 and thereby enable us to calculate the electrode- 
potential-dependent rate constants at various adatomic 
configurations as well as concentrations. Figure 01(c) 
demonstrates the rate constant of the reaction of ox¬ 
idation of ferrocyanide ions, k c , at graphene electrode 
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with potassium impurity as a function of electrode po¬ 
tential. At a potential of about —0.1 V vs. Ag/AgCl, 
increase of the cathodic reaction rate is observed. Con¬ 
trary to metallic electrodes, the increase of the k c is not 
monotonic. In the range of calculated potentials, the 
plot of k c vs. V has a hump. The local minimum ap¬ 
pears within the electrochemical window of water, i.e. 
within the range 0.6—0.8 V vs. Ag/AgCl. The monotonic 
increase is observed at the positive electrode potentials 
beyond the water window (V > 0.8 V vs. Ag/AgCl). Dif¬ 
ference between the k c vs. V plots calculated for the pure 
graphene and graphene with K adatoms is seen clearly. 
However, as well as the shape of DOS-curves close to 
the Dirac point [Fig. E[a)], the shape of fc c -pl°ts of im¬ 
pure graphene [Fig. |U[c)] is only negligibly affected by 
the impurity configuration. Generally, K-impurity slow 
down the reaction kinetics in the electrode potentials in 
the range of water window. At the pure graphene elec¬ 
trode, the maximum of the hump is located at poten¬ 
tial of 0.503 V vs. Ag/AgCl (k c = 0.194 a.u.), whereas 
at graphene with K adatoms in the concentration of 
3.125% the maximum is observed at lower potential of 
0.453 V vs. Ag/AgCl ( k c = 0.123 a.u.). The same applies 
to the position of the local minimum, which is down¬ 
shifted in graphene with K adatoms by 0.17 V. 

Contrary to the impurity configuration, the concen¬ 
tration of adatoms influences strongly the k c [Fig. [2(d)] 
similarly to the influence on the DOS [Fig. [2(b)]. With 
increasing adatomic concentration, the position of the 
local maximum shift towards lower potentials from the 
value of 0.503 V vs. Ag/AgCl (k c = 0.123 a.u.) for 
pure graphene to 0.393 V vs. Ag/AgCl (k c = 0.091 a.u.). 
Also the local minimum shifts down form the potential 
value of 0.843 V vs. Ag/AgCl to 0.393 V vs. Ag/AgCl 
to 0.573 V vs. Ag/AgCl for pure graphene and graphene 
with adatomic concentration of 12.5%, respectively. In 
a part of the k c plot in the range of higher potentials 
( V > 0.8 V vs. Ag/AgCl) an additional hump is appar¬ 
ent at 12.5% of adatoms. 

Assuming 0.2 V vs. Ag/AgCl as a standard electrode 
potential (when the rates of both cathodic and anodic re¬ 
actions are equal), in Table Hill we compare values of the 
standard rate constant, k s , for electron transfer between 
the Fe(CN)g redox couple and impure graphene at 
weakly long-range scattering potential in Fig. [3](a) and 
strongly short-range Gaussian potential V = Ue~ r 
with potential height U = —2 u and effective potential 
radius £ = 0.5a. While the adatomic configurations do 
not affect significantly the shape of the k c plot, apparent 
differences in the k s at different ranges of the scattering 
potential action are seen. When the long-range potential 
is used, the electron transfer is moderately suppressed 
by adatoms in random («31%), correlated (~33%), and 
ordered (~32%) configurations as compared to the elec¬ 
tron transfer of pure graphene. On the other hand, weak 
dependence between impurity concentration and the k s 
value is observed if the short-range potential is used; qua¬ 
drupling the adatomic content causes the decrease of k s 


Table III: Standard rate constants, k 3 , for electron trans¬ 
fer between graphene electrode and ferro-/ferricyanide redox 
couple for different ranges of the scattering-potential action, 
adatomic configurations and concentrations. 


Type of potential Configuration Stoichiometry 

k s [a.u.] 

Long-range 

Random 

Correlated 

Ordered 

1/32 (3.125%) 

0.0595 

0.0571 

0.0586 



1/32 (3.125%) 

0.0791 

Short-range 

Ordered 

1/16 (6.25%) 
1/8 (12.5%) 

0.0740 

0.0610 


Pure 

0 

0.0858 


by only «23%. Thus, the use of the long-range potential 
more strongly suppresses the electron transfer kinetics. 
It is worth noting that the best kinetics is observed at 
pure electrode. Our findings are not compatible with ex¬ 
perimental results obtained in Ref. |68]. Discrepancies are 
probably due to the hydrophobic properties of graphene. 

In contrast to the case of randomly-arranged adatoms, 
when steady diffusive regime is reached for a relatively 
short time [Fig. (5[a)], in case of their correlation and 
especially ordering, a quasi-ballistic regime is observed 
during a long time as it is shown in Figs. [5](b) and (c). 
This (quasi-ballistic) behavior of diffusivity, D(t), indi¬ 
cates a very low scattered electronic transport, at which 
maximal value of D(t) is substantially higher for corre¬ 
lated and much more for ordered adatoms as compared 
with their random distribution. If the diffusive regime 
is not completely reached, the semiclassical conductivity, 
<7, cannot be in principle defined. However, we extracted 
a for the case of ordered adatoms using the highest D(t) 
when quasi-ballistic behavior turns to a quasi-diffusive 
regime with an almost saturated diffusivity coefficient. 

Figure [6] represents calculated conductivity (er) as a 
function of electron (n e > 0) or hole (n e < 0) con¬ 
centration, a = cr(n e ), for different positions (viz. H 1 
B , and T) and distributions (viz. random, correlated, 
and ordered) of adatoms in graphene. For visual conve¬ 
nience, we arranged the same (nine) curves in two groups: 
Figs. [6][a)-(c) demonstrate how correlation and ordering 
affect the conductivity for each of H, B , and T adsorp¬ 
tion types, while Figs. E[d)-(f) exhibit how these three 
types of sites influence on the conductivity for each of 
random, correlated (with maximal correlation lengths as 
listed in Table CD, and ordered adatomic distributions. 
The conductivity exhibits linear or nonlinear (viz. sub- 
linear) electron-density dependencies. The linearity of 
a = a(n e ) takes place at randomly-distributed potassium 
adatoms and indicates dominance of the long-range con¬ 
tribution to the scattering potential, while sublinearity 
occurs at nonrandom (viz. correlated and ordered) posi¬ 
tions of K adatoms and is indicative of the dominance of 
short-range component of the scattering potential. This 
is in accordance with many previous studies (see, e.g., 
Ref. 86 and references therein) in which pronounced lin¬ 
earity and sublinearity of cr = cr(n e ) are observed for 








Figure 5: (Color online) Time-dependent diffusivity within the energy range E £ [—0.5u, 0.5u] for random (a), correlated (b), 
and ordered (c) potassium adatoms located on hollow sites. 
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Figure 6: (Color online) Conductivity vs. the electron density for uk = 3.125% of random, correlated, and ordered potassium 
adatoms occupying hollow ( H ), bridge ( B ), or top ( T ) adsorption sites. Curves in upper and lower figures are the same, but 
grouped in a different way to distinguish configuration effects induced by correlation or ordering from those caused by difference 
in type of adsorption sites: H, B, or T. 


long-range scattering potential (appropriate for screened 
charged impurities ionically bond to graphene) and short- 
range potential (appropriate for neutral covalently bond 
adatoms), respectively. These results illustrate manifes¬ 
tation of contrasting scattering mechanisms for different 
spatial distributions of metallic adatoms. 

One can see from Figs. [6[a)-(c) that conductivities 
for correlated (<7 cor ) and ordered (<T or d) adatoms are 
dozens of times enhanced as compared with case of 
randomly-distributed adatoms (<r rn d). These enhance¬ 
ments (<7 cor /<7rnd and <T or d/o'rnd) depends on electron den¬ 
sity and type of adsorption sites. It is easy to deter¬ 


mine from Figs. Ela)-(c) that the ratio a COI /a rn d ranges 
as 2<cr^ r (n e )/cr r ^ d (n e )<5, 2 < ag T (n e )/a? nd (n e ) < 6, 
and 3<CT^ r (n e )/a^ d (n e )<7; while cr md /<j m d ranges 
as 3 ^orf( n e)M® d W^8, 3<crf rd (n e )/cr r ^ d (n e )<9, 
and4<crj rd (n e )/cr 1 ^ d (n e ) < 15 (here, superscripts denote 
types of adsorption sites). As follows from Figs. (H[d)- 
(f), in a random adatomic state crJ nA < a^ nd < cr^ ld , while 
for correlated and ordered states, < a^ OI ~ a^ OI and 
CT ord ~ CT OTd ~ CT ord: respectively. That is why the highest 
increase of a due to correlation (<7 CO r/<7md) or ordering 
(oWd/cmd) takes place for the T-site bonding, followed 
by the B sites, and then the H sites. The increasing of 
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Figure 7: (Color online) Electron-density-dependent conductivity for different adsorption heights, h, of 3.125% of random (a), 
correlated (b), and ordered (c) K adatoms resided on hollow sites. 


cr due to adatomic correlation or ordering is expected in 
a varying degree for any constant-sign (V > 0 or V < 0) 
scattering potential, but this is not the case when the 
potential is sign-changing ( V 

If adatoms are randomly-positioned on the ( H , B , or 
T) adsorption sites, the conductivity is dependent on 
their type: a* d > a* d > aJ nA , particularly cr" d «2 f J 1 T nd 
[Fig-Eld)]. Here, the differences in a are caused by dif¬ 
ferent values of on-site potentials for these three types 
of adatomic positions although the same potential pro¬ 
file V = V{r) [Fig. [2[a)] is used for them. The stronger 
(weaker) on-site potential V) corresponds to the smaller 
(larger) distance r from the given graphene-lattice site 
i to the nearest adsorption site, which is more close 
(distant) for the H (T) type, followed by the B type, 
and then the T ( H ) type. If adatoms are correlated, 
the conductivity is dependent on whether they act as 
interstitial (H or B sites) or substitutional (T sites) 
atoms: a^ or « a^ OI > a^ or [Fig. ED(c)], which can be at¬ 
tributed to the values of maximal correlation lengths 
(Table [Til) defining correlation degree for H 1 B and T 
sites, rff =r£ > rJ . Finally, if adatoms form or- 

' '-'max '-'max '-'max " ' 

dered superstructures (superlattices) with equal periods 
(Fig. [3]), the conductivity is practically independent on 
the adsorption type (especially for not very high charge 
carrier densities): cr£ d w of rd « crf rd [Fig. G3T)] . 

In our model, increase (or decrease) of adatomic eleva¬ 
tion over the graphene surface results to more weak (or 
strong) scattering-potential amplitude, i.e. physically it 
means more weak (or strong) regime of electron scatter¬ 
ing on charged impurity adatoms. Although the values of 
adsorption height, h, reported in the literature for potas¬ 
sium do not disagree as much as for the adsorption energy 
(see TablelU), for the model and calculation completeness, 
we range h in a wide interval (up to h = 3.6 A) includ¬ 
ing an exotic case of h = 0, when impurity atoms act 
as strictly interstitial ones. Calculated curves represent¬ 
ing the charge-carrier-density-dependent conductivity for 
(random, correlated, and ordered) adatoms resided on 
(the most favorable for potassium) hollow sites and ele¬ 
vated on different h are shown in Fig. [TJ (Here, we do not 


consider the cases of less favorable for potassium bridge 
and top sites since it leads to qualitatively the same re¬ 
sults.) As follows from Figs. [7[a) and (b), at least for 
hole densities (— n e > 0), two (three) time increased or 
decreased h for randomly- or correlatively-distributed K- 
adatoms results to approximately two (three) time en¬ 
hanced or reduced a, respectively. Thus the conductiv¬ 
ity approximately linearly scales with adsorption height 
of random or correlated adatoms, cr{h) oc h. However, 
for ordered potassium adatoms, the cr remains practi¬ 
cally unchanged with varying of h in the realistic range 
of adsorption heights (see Table [I]) and even in all range 
at issue (0 < h ^ 3.6 A) for hole densities [Fig. [TJ^c)] . 
We attribute this to the dominance of short-range scat¬ 
tered in case of their ordered state as it was mentioned 
above. Indeed, the Gaussian fitting for the scattering 
potential in Fig. [2[a) yields the effective potential radius 
£ = 2.21a, which is commensurable with quantities of ad¬ 
sorption heights h at issue (and even less than h = 3.6 A). 

In conclusion of this section, note that our numeri¬ 
cal calculations of conductivity in Figs. [6] and [T] agree 
with experimentally observed features of cr = a(n e ) in 
potassium-doped graphene^^ (i) asymmetry in the con¬ 
ductivity for electrons versus holes (which, however, can 
be weakened and even totally suppressed due to the spa¬ 
tial correlation or especially ordering of adatoms as well 
as increasing of their adsorption height), (ii) shifting of 
minimum conductivity at a charge neutrality point to 
more negative gate voltage, (iii) linearity or sublinear- 
ity of conductivity at lower or higher gate voltage, re¬ 
spectively, and (iv) increase in conductivity due to cor¬ 
relation in the positions of adatoms that was also sus¬ 
tained theoretically^ within the standard semiclassical 
Boltzmann approach in the Born approximation. A sig¬ 
nificant sublinear behavior of electron-density-dependent 
conductivity and its saturation for very high densities at 
the spatial correlations among the charged impurity loca¬ 
tions in contrast to the strictly linear-in-density graphene 
conductivity for uncorrelated random charged impurity 
scattering [Figs. [6[a)~(e) and Efa)-(b)] is also in agree¬ 
ment with theoretical findings in Refs. HH andf99l 
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IV. CONCLUSIONS 

By employing numerical calculations, we systemati¬ 
cally studied the effects of different (random, correlated, 
and ordered) spatial configurations of potassium adatoms 
onto high-symmetry [hollow- ( H ), bridge- (T), and top- 
type ( T)} adsorption sites with various elevations over 
the graphene sheet on its electron transport and elec¬ 
trochemical properties to ascertain correlation between 
them. We conclude as follows. 

(i) The charge carrier density dependence of the con¬ 
ductivity is indicative of dominance of long-range scat¬ 
tering centers for their random spatial distribution, while 
short-range scatterers dominate for their correlated and 
ordered states. This demonstrates manifestation of con¬ 
trasting scattering mechanisms for different spatial dis¬ 
tributions of metallic adatoms. 

(ii) A band gap may be opened only if ordered adatoms 
act as substitutional atoms ( i.e. reside on T-type sites) 
due to the breaking of graphene lattice point symmetry, 
while there is no band gap opening for adatoms acting 
as interstitial atoms (i.e. occupying H- or 5-type sites). 

(iii) If adatoms are randomly-positioned on the H , B- 
or T sites, the conductivity is dependent on their type: 
a rnd > °rnd > °md- For spatially-correlated adatoms, the 
conductivity is dependent on whether they act as in¬ 
terstitial or substitutional atoms: a^ OT « a^ OI > a^ or . If 
adatoms form ordered superstructures (superlattices) 
with equal periods, the conductivity is practically inde¬ 
pendent on the adsorption type (especially for low elec¬ 
tron densities): <j^ d « crf rd « crf rd . 

(iv) Depending on electron density and type of ad¬ 
sorption sites, the conductivity for correlated and or¬ 
dered K adatoms is found to be enhanced in dozens of 
times as compared to the cases of their random positions. 
The correlation and ordering effects manifest stronger for 
adatoms acting as substitutional atoms and weaker for 
those acting as interstitial atoms. 

(v) The electron-hole asymmetry in the conductivity 


for randomly-positioned adatoms weakens and even may 
be totally suppressed for correlated and especially or¬ 
dered ones as well as for increased of their adsorption 
height. 

(vi) The conductivity dependence with adsorption 
height of random or correlated adatoms scales approx¬ 
imately as a(h) oc h. However, for ordered adatoms, a 
remains practically unchanged with varying of h within 
its realistic range. 

(vii) Only slight suppress of electron transfer kinet¬ 
ics in electrolyte at K-doped graphene electrode is re¬ 
vealed. Strong correlation between the band gap in 
graphene and the shape of the electrode-potential depen¬ 
dence of electrochemical rate constant is seen, when the 
strongly short-range scattering potential during the elec¬ 
tron transport in graphene is used. At the same time, the 
influence of this potential on the suppress of the standard 
electrochemical rate constant is much weaker as com¬ 
pared to the case of the long-range electron-scattering 
potential in graphene. Comparison of the electron trans¬ 
fer calculations to experiment shows that the hydropho- 
bicity of graphene is a key factor, which suppresses the 
kinetics of heterogeneous electron transfer in electrolyte 
at graphene electrode. 
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